Active nematics are intrinsically phase-separated 
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Two-dimensional nonequilibrium nematic steady states, as found in agitated granular-rod mono- 
layers or films of orientable amoeboid cells, were predicted [Europbys. Lett. 62 (2003) f 96] to have 
giant number fluctuations, with standard deviation proportional to the mean. We show numerically 
that the steady state of such systems is macro scopically phase-separated, yet dominated by fluctu- 
ations, as in the Das-Barma model [PRL 85 (2000) f602]. We suggest experimental tests of our 
findings in granular and living-cell systems. 

PACS numbers: 05.70.Ln Nonequilibrium and irreversible thermodynamics. 87. 18. Ed Aggregation and 
other collective behavior of motile cells. 45.70.-n Granular systems. 



The ordering or "flocking" Q, 0, Q of self-propelled 
particles obeys laws strikingly different from those 
governing thermal equilibrium systems of the same 
spatial symmetry. Even in two dimensions, the ve- 
locities of particles in such flocks show true long- 
range order 0, Q] , despite the spontaneous breaking 
of continuous rotational invariance. Density fluctua- 
tions in the ordered phase are anomalously large Q , 
and the onset of the ordered phase is discontinuous 
. The ultimate origin of these nonequilibrium phe- 
nomena is that the order parameter is not simply an 
orientation but a macroscopic velocity. It is thus in- 
triguing that even the nematic phase of a collection 
of self-driven particles, which is apolar and hence has 
zero macroscopic velocity, shows 0, giant number 
fluctuations |7j, as a result of the manner in which 
orientational fluctuations drive mass currents. This 
Letter takes a closer look at these fluctuations and 
shows that they offer a physical realisation of the 
remarkable nonequilibrium phenomenon known as 
fluctuation- dominated phase separation 0, hitherto 
a theoretical curiosity. 



serves as the motor coordinate). An active nematic 
is a collection of such particles with axes on average 
spontaneously aligned in a direction n, with invari- 
ance under n — > n. We know of two realisations of 
active nematics: collections of living amoeboid cells 
[T(ij and granular-rod monolayers 0, . 
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FIG. 1: Number standard deviation An scales roughly 
as the mean N, for system sizes L — 128, 256 



Before presenting our results, we make precise the 
term active nematic. An active particle extracts en- 
ergy from sources in the ambient medium or an in- 
ternal fuel tank, dissipates it by the cyclical motion 
of an internal "motor" coordinate, and moves as a 
consequence. For the anisotropic particles that con- 
cern us here, the direction of motion is determined 
predominantly by the orientation. Our definition 
encompasses self-propelled organisms, living cells, 
molecular motors, and macroscopic rods on a ver- 
tically vibrated substrate (where the tilt of the rod 
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We study active nematics in a simple numerical 
model described in detail below. Our results confirm 
(see Fig. the giant number fluctuations (standard 
deviation oc mean) [(| predicted by the linearised 
analysis of but are far richer: (i) A statisti- 
cally uniform initial distribution of particles, on a 
well-ordered nematic background, undergoes a del- 
icate "fluctuation-dominated" @ phase separation, 
where the system explores many statistically sim- 
ilar segregated configurations, (ii) The equal-time 
two-point density correlator C(r), Fig. |5] shows a 
collapse when plotted as a function of r/L(t), where 
L{t) is the location of the first zero-crossing of C(r), 
with a cusp at small r/L(t) signalling a departure 
from Porod's Law, i.e., the absence of sharp inter- 
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FIG. 2: Coarsening length L(t) saturates to a value pro- 
portional to system size. Inset: L(t) consistent with t 1 ' 2 
despite conservation law. 



faces, (iii) We confirm that the large density in- 
homogeneities are indeed best thought of as phase 
separation, by showing (1) that the saturation value 
of L(t) is proportional to the linear size of the sys- 
tem and (2) that the phase-separation order param- 
eter - the time-averaged first spatial Fourier compo- 
nent of the particle density - approaches a nonzero 
value in the limit of large system size (Fig. 0) . (iv) 
L(t) grows clearly faster than the i 1 / 3 for normal 
conscrvcd-ordcr-parameter coarsening (see inset to 
Fig. 0), and is consistent with i 1 / 2 as expected by 
analogy to . Below we show how these results were 
obtained, discuss them in detail, explain the analogy 
to the work of || , and suggest experimental tests for 
these striking phenomena. 



tors but equal-time quantities like Sq or An as wefl 
must be inferred from equations of motion for n and 
c. As shown in Q , the equation of motion for n is 
qualitatively the same as for equilibrium ncmatics. 
The feature [f| that distinguishes active nematics 
crucially from their equilibrium counterparts is that 
the current j in the continuity equation d t c = — V • j 
for the density has a contribution oc V • c(nn) [l4j . 
This term, which is ruled out at thermal equilib- 
rium, has a simple, physically appealing origin: spa- 
tial variation in the director field n defines a curve; 
the normal to this curve defines a local vectorial 
asymmetry; for a driven system, such an asymme- 
try implies a current 0|. If n = (cos sin 9) then 
inhomogencitics in 9 give a curvature-induced cur- 
rent j = (jx,jz) oc (d6/dz,d9/dx) (see Fig. EJ, 
analogous to Q where particles slide with velocity 
oc V/i on a fluctuating interface with height field 
h, except that our current is not a gradient, the 
Since nematic order is a spontaneous breaking of 
rotation-invariance, large fluctuations in 9 at long 
wavelengths are expected to be present in abun- 
dance, and to decay slowly, in any nematic, equi- 
librium or otherwise. We showed in the previous 
paragraph that the effect of these broken-symmetry 
modes on the density field was benign in an equi- 
librium nematic. In an active nematic, however, 
the same orientational fluctuations, because of the 
curvature-induced current we just mentioned, will 
affect the density fluctuations substantially. A lin- 
earised small-fluctuations analysis [j| showed that 
they to lead to giant fluctuations in the number of 
particles: Ajv/V^/V" oc N 1 ^ in d dimensions, i.e., 
A N oc TV for d = 2. 



We begin with some background information. An 
apolar, uniaxial, compressible nematic liquid crystal 
is described by director and number-density fields 
n(r) and c(r), with fluctuations <5n(r) and 5c(r) 
about their uniform mean values no and cq. Let 
us review first what happens at thermal equilibrium. 
The system is then governed by an extended Frank 
[3 free-energy F[n,c] = (l/2)/d d r[K(Vn) 2 + 
A(Sc) 2 /c + dn ■ VcV • n + C 2 n x Vc • V x n], 
where K is an elastic tensor, A the comprcssional 
modulus at constant orientation, and C1.2 couple 
orientation and density in the simplest symmetry- 
allowed fashion. Equipartition applied to F implies 
that the static structure factor £ q = J d d rexp(— iq- 
r)(<5c(0)<5c(r))/crj is finite for q — ► 0; i.e., the mean 
TV and standard deviation An of the number TV of 
particles obey An oc V / TV at equilibrium even when 
Ci, c 2 ^0. 

An active nematic is a steady state away from ther- 
mal equilibrium, in which not only dynamic correla- 



Such large fluctuations prompt the suspicion that an 
analysis beyond Gaussian fluctuations would reveal 
that the system is in fact phase-separated, as in ||. 
There are two issues here: (i) whether the nonequi- 
librium coupling mentioned above inevitably arises 
in an active nematic; and (ii) whether it leads to 
phase separation. Ref. @ effectively answers the 
first question in the affirmative; we focus on the sec- 
ond. 

We find it convenient to separate the density and 
orientation degrees of freedom, and employ a dis- 
crete model of lattice-gas particles coupled to an an- 
gle field, incorporating explicitly the the nonequi- 
librium curvature-induced current j oc (d z 9, d x 9) 
mentioned above, via a suitable choice of particle- 
hopping rates. We consider a two-dimensional lat- 
tice with angles 9i e [0, tt] and noninteracting lattice- 
gas occupancy variables rii = 0, 1 at each site i. 
The angles evolve by Metropolis Monte Carlo up- 
dates governed by the Lebwohl-Lasher 0] hamilto- 
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FIG. 3: Director variation from the bottom to the top of 
the picture gives rise to a current in a transverse direction 



nian H = —K^2 <i - > cos2(9 i — 8j) yielding a ne- 
matic phase at low temperature. Particle motion 
is nonequilibrium: Hops of a particle at site i to a 
nearest-neighbour site in the ±x direction are at- 
tempted with probability 1/4 ± a(6x — 6%), and in 
the ±z direction with probability 1/4 ± 0.(83 — 6*4), 
where Oi are the angles at sites i = 1 to 4 as in 
Fig. |21 and a encodes the strength of the ac- 
tive curvature-current coupling of 0. In the 




FIG. 4: Anisotropic clustering of particles, 20000 time 
steps, system size 100, 50 % occupancy 

effect of the concentration field on the dynamics of 
the angle field is shown to play an insignificant role 
in the giant number fluctuations. Accordingly, we 
neglect it here, so that the particles are advected 
passively by an autonomous angle field. This simpli- 
fying approximation should not make a qualitative 
difference deep in the nematic phase. Accordingly, 
we work with K = 10.0, and let the angle field equi- 
librate until a well-ordered nematic is formed. An 
initial statistically uniform distribution of particles, 
(mean occupancy —50 %) clusters and coarsens (Fig. 
BJ, most strongly at ±45° to the mean direction of 
nematic ordering. The anisotropic two-point den- 
sity correlator C(r,t) shows a scaling collapse for a 
given direction for all t, if plotted as a function of 
r/L(t) for a coarsening length L(t), defined by the 
value of r at the first zero-crossing, whose scaling 
is consistent with t 1 / 2 . (see Fig. and Fig. |2J). 
The value at which L(i) saturates is proportional 
to the system size L (see Fig. which makes a 
strong case for true phase separation. The value of 



the saturation length is numerically small comparer} 
to i, probably because of the poorly-defined clus- 
ters (Fig. 0J) of fluctuation-dominated phase separa- 
tion. That we work with hard-core particles, and on 
timescales on which the macroscopic variation of the 
mean nematic orientation is very small, also proba- 
bly contributes. The exponent of 1/2 is because the 
particles aggregate not by diffusion plus short-range 
capture, but rather, by analogy with Q (see also 
|la | ) , because the broken-symmetry mode of the ne- 
matic order sweeps the particles over large distances 
via curvature-induced drift. A nematic fluctuation 
on a scale £ collects particles in a time of order i. For 
two such domains to coalesce requires the nematic 
director field on that scale to turn over, which is a 
time ~ t z where z = 2 is the dynamic exponent of 
transverse fluctuations of the nematic director. 
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FIG. 5: Plots of the equal-time density correlator C(r,t) 
at different times t collapse onto a single curve when r is 
scaled by the zero-crossing coarsening length L(t); data 
shown for direction transverse to nematic ordering Inset: 
C(r,t) vs. r. 



At long times a steady state is reached, and 
C(r/L(t — > 00)) shows a cusp at small argument 
(C(x) ex x a ,a ~ 0.33), which can be seen in Fig. 
as well, signalling the absence of sharp interfaces 
between regions rich and poor in particles, and a 
power-law distribution of cluster sizes. For steady 
state in the largest system, we also measured the 
standard deviation Ajy in the number of particles 
in an observation containing N particles on average. 
The plot of An vs N, Fig. ^ shows precisely the 
linear dependence predicted by || . Faced with these 
results, we ask: is this phase separation or a single 
phase with large fluctuations? This is answered by 
measuring the magnitude of the time-averaged low- 
est spatial Fourier-component Q(l, 1) of the density, 
shown in Fig. |SJ Although the data are not con- 
clusive, the flattening of the semilog plot as a func- 
tion of system size rules out an exponential decay 
to zero. Together with the proportionality of the 



coarsening length to the system size, and the nature 
of the mechanism, this suggests strongly that active 
ncmatics offer the most natural physical realisation 
of macroscopic fluctuation-dominated phase separa- 
tion ||. As in 0, we find that the time-series of 
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FIG. 6: Time-averaged fourier component of the density 
profile vs. wave vector k for wave vector direction (1,1), 
for system sizes L = 16, 32, 64, 128; flattening in semilog 
plot at largest size suggests nonzero value for L — > oo. 
Inset: Time series of Q for (1,1) direction. 

Q(l, 1) shows enormous fluctuations, (Fig. El inset), 
with "crashes" during which other nearby fourier- 
components gain weight. Thus, as in the system 
lurches from one macroscopically phase-separated 



configuration to another, spending very little time irl 
non-phase-separated states. Lastly, the velocity au- 
tocorrelation of tagged particles at low (15 %) con- 
centration agrees qualitatively with the l/t tail (plot 
not shown) predicted by [f| , over the range in which 
a given particle moves unimpeded by others. 

What experiments can test these results? The best 
would be agitated layers of granular rods, for which 
nematic phases have been reported Although 
many features of these systems can be rationalized 
in terms of equilibrium hard-rod theories 01 , some 
properties such as global circulation and swirls fill 
|l2J are clearly very noncquilibrium. These systems 
as well as the living melanocyte nematic of Grulcr 
et al. [Tol | remain the most promising candidates for 
experimental tests of the rich range of results made 
here and in p|. The confirmation of giant number 
fluctuations in the numerical experiments of is 
encouraging in this regard. 
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